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Abstract 

The hard-core model has received much attention in the past couple of decades as a 
lattice gas model with hard constraints in statistical physics, a multicast model of calls 
in communication networks, and as a weighted independent set problem in combinatorics, 
probability and theoretical computer science. 

In this model, each independent set / in a graph G is weighted proportionally to A'^', 
for a positive real parameter A. For large A, computing the partition function (namely, the 
normalizing constant which makes the weighting a probability distribution on a finite graph) 
on graphs of maximum degree A > 3, is a well known computationally challenging problem. 
More concretely, let Ac(Ta) denote the critical value for the so-called uniqueness threshold 
of the hard-core model on the infinite A-regular tree; recent breakthrough results of Dror 
Weitz (2006) and Allan Sly (2010) have identified Ac(Ta) as a threshold where the hardness 
of estimating the above partition function undergoes a computational transition. 

We focus on the well-studied particular case of the square lattice Z^, and provide a new 
lower bound for the uniqueness threshold, in particular taking it well above Ac(T4). Our 
technique refines and builds on the tree of self- avoiding walks approach of Weitz, resulting 
in a new technical sufficient criterion (of wider applicability) for establishing strong spatial 
mixing (and hence uniqueness) for the hard-core model. Our new criterion achieves better 
bounds on strong spatial mixing when the graph has extra structure, improving upon what 
can be achieved by just using the maximum degree. Applying our technique to I? we 
prove that strong spatial mixing holds for all A < 2.3882, improving upon the work of 
Weitz that held for A < 27/16 = 1.6875. Our results imply a fully-polynomial deterministic 
approximation algorithm for estimating the partition function, as well as rapid mixing of the 
associated Glauber dynamics to sample from the hard-core distribution. 
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1 Introduction 



In this paper we study phase transitions for samphng weighted independent sets (weighted by 
an activity A > 0) of the 2-dimensional integer lattice 1? . In statistical physics terminology, we 
study the hard-core lattice gas model ([ , ]), which is a simple model of a gas whose particles 
have non-negligible size (thus preventing them from occupying neighboring sites), with activity 
A G M+ corresponding to the so-called fugacity of the gas. More formally, for a finite graph 
G = iy^E)^ let = fi(G) denote the set of independent sets of G. Given an independent set 
(J G ri, its weight is defined as w{g^ — A'^' and v \^ said to be occupied under g \i v ^ a. 
The associated Gibbs (or Boltzmann) distribution — fiQ^x is defined on Q as fi(a) = w(a)/Z, 
where Z = Z(G, A) = ^r]en ^(v) is commonly referred to as the partition function. 

Recall that Valiant [ ] showed that exactly computing the number of independent sets is 
#P-complete, even when restricted to 3-regular graphs (see Greenhill [ ]). Hence, we focus 
our attention on approximation algorithms for estimating the number, or more generally, the 
partition function. It is well known [ ] that the problem of approximating the partition func- 
tion Z and that of sampling from a distribution that is close to the Gibbs distribution /x, are 
polynomial-time reducible to each other (see also [ ]). 

The fundamental notion of a phase transition for a statistical mechanics model on an infi- 
nite graph addresses the critical point at which the model starts to exhibit a certain long-range 
dependence, as a system parameter is varied. In particular, the so-called critical inverse tem- 
perature /3c for the Ising or the Potts model, and the critical activity Ac for the hard-core lattice 
gas model, are prime examples where the system undergoes a transition from uniqueness to 
multiplicity of the infinite- volume Gibbs measures. 

Phase transition in the hard-core model is also intimately related to the computational com- 
plexity of estimating the partition function Z. Recently, a remarkable connection was estab- 
lished between the computational complexity of approximating the partition function for graphs 
of maximum degree A and the phase transition Ac(Ta) for the infinite regular tree Ta of degree 
A. On the positive side, Weitz [34] showed a deterministic fully-polynomial time approximation 
algorithm (FPAS) for approximating the partition function for any graph with maximum degree 
A, when A < Ac(Ta) and A is constant. On the other side. Sly [ ] recently showed that for 
every A > 3, it is NP-hard (unless NP=RP) to approximate the partition function for graphs of 
maximum degree A, when Ac(Ta) < A < Ac(Ta) + eA, for some function e/\ > 0. More recently, 
Galanis et al. [12] improved the range of A in Sly's inapproximability result, extending it to all 
A > Ac(Ta) for the cases A = 3 and A > 6. 

1.1 Prior history and current work 

Our work builds upon Weitz's work to get improved results for specific graphs of interest. We 
focus our attention on what is arguably the simplest, not yet well-understood, case of interest 
namely the square grid, or the 2-dimensional integer lattice Z^. Empirical evidence suggests 
that the critical point Xc{Z'^) ^ 3.796 [13, 3, 26], but rigorous results are significantly far from 
this conjectured point. The possibility of there being multiple such Ac is not ruled out, although 
no one beheves that this is the case. 

From below, van den Berg and Steif [ ] used a disagreement percolation argument to prove 
Ac(Z^) > where Pc is the critical probability for site percolation on Z^. Applying the best 
known lower bound on pc > 0.556 for 7? by van den Berg and Ermakov [5] implies Ac(Z^) > 
1.252.... Prior to that work, an alternative approach aimed at establishing the Dobrushin- 
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Shlosman criterion [10], yielded, via computer-assisted proofs, Ac(Z^) > 1.185 by Radulescu and 
Styer [ ], and Ac(Z2) > 1.508 by Radulescu [27]. 

These results were improved upon by Weitz [ ] who showed that Ac(Z^) > Ac(T4) = 27/16 = 
1.6875, where Ta is the infinite, complete, regular tree of degree A. For the upper bound, a 
classical Peierls' type argument implies Ac(Z^) = 0(1) [ ]. (A related result of Randall [ ] 
showing slow mixing of the Glauber dynamics for A > 8.066 gives hope for a better upper bound 
on Ac(Z2).) The regular tree Ta is one of the only examples (that we know of) where the critical 
point is known exactly, and in this case, Kelly [ ] showed that Ac(Ta) = (A - 1)^-V(A - 2)^. 

In this work we present a new general approach which, for the case of the hard-core model 
on Z^, improves the lower bound to Ac(Z2) > 2.3882. There are various algorithmic implications 
for finite subgraphs of the 1? when A < 2.3882. Our results imply that Weitz 's deterministic 
FPAS is also valid on subgraphs of 1? for the same range of A. Thanks to the existing literature 
on general spin systems (['^'^, 23, 8, 11]), our results also imply that the Glauber dynamics has 
O(nlogn) mixing time for any finite subregion G — (y, E) of 1? when A < 2.3882, where n— \y\. 
Recall that the Glauber dynamics is a simple Markov chain that updates the configuration at 
a randomly chosen vertex in each step, see [ ] for an introduction to the Glauber dynamics. 
The stationary distribution of this chain is the Gibbs distribution. Hence, it is of interest as 
an algorithmic technique to randomly sample from the Gibbs distribution, and also as a model 
of how physical systems reach equilibrium. The mixing time is the number of steps (from the 
worst initial configuration) until the distribution is guaranteed to be within variation distance 
< 1/4 of the stationary distribution. 

As in Weitz's work, our approach can be used for other 2-spin systems, such as the Ising 
model. This is discussed in Section 6. As will be evident from the following high-level idea of 
our approach, it can be applied to other graphs of interest. Our work also provides an arguably 
simpler way to derive the main technical result of Weitz showing that any graph with maximum 
degree A has strong spatial mixing (SSM) when A < Ac(Ta). 

To underline the difficulty in estimating bounds on Ac, we remark that the existence of a 
(unique) critical activity Ac remains conjectural and an open problem for Z^, for d > 2. In 
contrast, for the Ising model, the critical inverse temperature /3c(Z^) has been known since 1944 
[24]; interestingly, the corresponding critical point for the g'-state Potts model (for > 2) has 
only recently been established (by Beffara and Duminil-Copin [ ]) to be /3c(^) = log(l + y^), 
settling a long-standing open problem. The lack of monotonicity in A in the hard-core model 
poses a serious challenge in establishing such a sharp result for this model. In fact, Brightwell 
et al. [ ] showed that in general such a monotonicity need not hold, by providing an example 
with a non-regular tree. 

2 Technical Preliminaries and Proof Approach 

Before presenting our approach, it is useful to review briefly the uniqueness/non- uniqueness 
phase transition, and introduce associated notions of decay of spatial correlation, known as 
weak and strong spatial mixing properties. Much of the below discussion is simplified for the 
case of the hard-core model on Z^, wherein one utilizes certain induced monotonicity (given by 
the bipartite property) in the model and the amenability of the graph. 
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2.1 Uniqueness, Weak and Strong Spatial Mixing 

Let Bl denote the finite grapli corresponding to a box of side-length 2L + 1 centered around 
the origin in Z^. Thus, Bl = {V, E), where V = (i, j) G I? : -2L - 1 < i, j < 2L + 1 with edges 
between pairs of vertices at Li distance (or Manhattan distance) equal to one. Since this is a 
bipartite graph, we may fix one such partition V — even U odd - for example, it is standard to 
consider the set of vertices at an even distance from the origin as the even set. The boundary 
oi Bl are those vertices v — {vi^V2) where \vi\ = 2L + 1 ioi i — 1 or z = 2. The hard-core 
model on bipartite graphs is a monotone system (e.g., see [ ]), which for the current discussion 
implies that we only have to consider two assignments to the boundary: all even vertices or all 
odd vertices on the boundary are occupied. Let a^^f^ (<^L^r ) denote the marginal probability 
that the origin r is unoccupied given the even (odd, respectively) boundary. Then to establish 
uniqueness of the Gibbs measures, we need that: 

lim |«iX-<'l-0- 

We are interested in the critical point Ac for the transition between uniqueness and non- 
uniqueness. A standard way to establish uniqueness is by proving one of the spatial mixing 
properties introduced next. 

Let G = (y, E) be a (finite) graph. For 5 C ^ , a configuration p on 5 specifies a subset of 
S as occupied and the remainder as unoccupied. Let /x^ = /x^ denote the Gibbs distribution 
conditional on configuration p to S. For v ^ let ay = aQ ^ denote the marginal probability 
that V is unoccupied in fi^. 

The first spatial mixing property is Weak Spatial Mixing (WSM). Here we consider a pair of 
boundary configurations on a subset S and consider the "influence" on the marginal probability 
that a vertex v is unoccupied. WSM says that the influence on v decays exponentially in the 
distance of S from v. 

Definition 1 (Weak Spatial Mixing). For the hard-core model at activity for finite graph 
G = (V^E), WSM holds with rate 7 G (0, 1) if for every v ^ V, every S C V, and every two 
configurations p, r] on 

where dist{v, S) is the graph distance (i.e., length of the shortest path) between v and (the nearest 
point in) the subset S. 

The second property of interest is Strong Spatial Mixing (SSM). The intuition is that if a 
pair of boundary configurations on a subset S agree at some vertices in S then those vertices 
"encourage" v to agree. Therefore, SSM says that the influence on v decays exponentially in 
the distance of v from the subset of vertices where the pair of configurations differ. 

Definition 2 (Strong Spatial Mixing). For the hard-core model at activity X, for finite graph 
G = (V, E), SSM holds with rate 7 G (0, 1) if for every v , every S dV , every S' C S, and 
every two configurations p^rj on S where p(S \ S') = r](S \ S'), 

Note that since dist(i;, T) < dist(i;, T\S), SSM implies WSM for the same rate. Moreover, it 
is a standard fact that such an exponential decay in finite boxes (say), in Z^, implies uniqueness 
of the corresponding infinite volume Gibbs measure on Z^, see Georgii [ ] for an introduction 
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to the theory of infinite-volume Gibbs measures. We can speciahze the above notions of WSM 
and SSM to a particular vertex in which case we say that WSM or SSM holds at v. If the 
graph is a rooted tree, we will always assume that the notions of WSM and SSM are considered 
at the root. 

For the hard-core model on a graph G = {V^E), for a subset of vertices S and a fixed 
configuration p on 5, it is equivalent to consider the subgraph which we obtain for each 

V ^ S that is fixed to be unoccupied we remove v from G, and for each v ^ S that is fixed to 
be occupied we remove v and its neighbors N{v) from G. In this way we obtain the following 
observation which will be useful for proving SSM holds. 

Observation 1. For a graph G = {V, E) and v eV, SSM holds in G at vertex v iff WSM holds 
for all subgraphs G' (of G) at vertex v. To he precise, by subgraphs we mean graphs obtained by 
considering all subgraphs of G and taking the component containing v. 

2.2 Self- Avoiding Walk Tree Representation 

Since our work builds on that of Weitz's, we first describe the self-avoiding walk (SAW) tree 
representation introduced in [34]. Given G = {V, E), we first fix an arbitrary ordering >w on the 
neighbors of each vertex w in G. For each v ^V^ the tree Tsaw(G,i;) is constructed as follows. 
Consider the tree T of self- avoiding walks originating from additionally including the vertices 
closing a cycle as leaves of the tree. We then fix such leaves of T to be occupied or unoccupied 
in the following manner. If a leaf vertex closes a cycle in G, say w ^ vi ^ . . . V£ ^ then if 
'^1 >w "^i we fix this leaf to be unoccupied, otherwise if vi V£ we fix the leaf to be occupied. 
Note, if the leaf is fixed to be unoccupied we simply remove that vertex from the tree. If the leaf 
is fixed to be occupied, we remove that leaf and all of its neighbors, i.e. we remove the parent 
of that leaf from the tree. The resulting tree is denoted as Tgaw = 7saw(G, i;). See Figure 1 for 
an illustration of Tgaw ioi a particular example. 

Weitz [ ] proves the following theorem for the hard-core model, which shows that the 
marginal distribution at the root in rsaw(G, i;) is identical to the marginal distribution for v in 
G. For a graph G = {V,E), a subset S CV and configuration p on 5, for Tsaw — 2^saw(G, v^, let 
p in Tgaw denote the configuration on S in Tgaw where for w ^ S every occurrence of w in Tgaw 
is assigned according to p. 

Theorem 1 (SAW Tree Representation, Theorem 3.1 in [ ]). For any graph G = (V,E), 

V ^ V , X > 0, and configuration p on S dV , for T — Tsa,w{G^ v) the following holds: 

p p 

Note, the tree Tgaw(G, i;) preserves the distance of vertices from v in G, which implies the 
following corollary. 

Corollary 2. // SSM holds with rate 7 for Tgaw(G, i;) for all v, then SSM holds for G with 
rate 7. 

The reverse implication of Corollary 2 does not hold since there are configurations on S in 
Tgaw which are not necessarily realizable in G. Observe that if G has maximum degree A, any 
SAW tree of G is a subtree of the regular tree of degree A. 



5 




: Occupied 
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Fi gure 1: Example of self- avoiding walk tree Tsaw 

The above tree describes Tsaw(G,a) with 
occupied and unoccupied leaves, while the below one is the same tree after removing those 
assigned leaves. At each vertex, we consider the ordering N > E > S > W ofits neighbors where 
5, W represent the neighbors in the North, East, South, West directions, respectively. 



2.3 Our Proof Approach 

In summary, Weitz [ ] first shows (via Theorem 1) that to prove SSM holds on a graph G = 
(V, E), it suffices to prove SSM holds on the trees Tqq^-^ {G,v), for ah v eV. Weitz then proves 
that the regular tree T/\ "dominates" every tree of maximum degree A in the sense that, for all 
trees of maximum degree A, SSM holds when A < Ac(Ta). We refine this second part of Weitz's 
approach. In particular, for graphs with extra structure, such as G = Z^, we bound Tq^.^{I?) 
by a tree T* that is much closer to it than the regular tree Ta- We then establish a criterion 
that achieves better bounds on SSM for trees when the trees have extra structure. 

The tree T* will be constructed in a regular manner so that we can prove properties about 
it - the construction of T* is governed by a (progeny) txt matrix M ^ whose rows correspond to 
t types of vertices, with the entry Mij specifying the number of children of type j that a vertex 
of type i begets. We will then show a sufficient condition using entries of M which implies 
that SSM holds for T* and for any subgraph of T*, including rsaw(^^)- The construction of 
T* is reminiscent of the strategy employed in [ , ] to upper bound the connectivity constant 
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of several lattice graphs, including 1? . The derivation of our sufficient condition has some 
inspiration from belief propagation algorithms. 

As a byproduct of our proof that our new criterion implies SSM for T*, we get a new (and 
simpler) proof of the second part of Weitz's approach, namely, that for all trees of maximum 
degree A, SSM holds when A < Ac(Ta). 



3 Branching Matrices and Strong Spatial Mixing 

As alluded to above, we will utilize more structural properties of self-avoiding walk trees. To 
this end, we consider families of trees which can be recursively generated by certain rules; we 
then show that such a general family is also analytically tractable. 



3.1 Definition of Branching Matrices 

We say that the matrix M is a, txt branching matrix if every entry Mij is a non-negative integer. 
We say the maximum degree of M is A = A{M) — maxi<-^<^ ^-[^^^-^^ M^j, the maximum row 
sum. Given a branching matrix M, we define the following family of graphs. In essence, it 
includes a graph G if the self-avoiding walk trees of G can be generated by M. 

Definition 3 (Branching Family). Given atxt branching matrix M , T<m includes trees which 
can be generated under the following restrictions: 

o Each vertex in tree T E T<m has its type i E {1, . . . , t}. 

o Each vertex of type i has at most Mij children of type j. 

In addition, we use the notation G — (V, E) E J^<m if Tqq^^^{G^ v) E T<m for all v ^V. 

For example, the family J^<m with M = [A] includes the family of trees with maximum 
branching A. On the other hand, J^<m with M = ( ^ ^ ) describes the family of graphs 



,0 A 

of maximum degree A + 1, by assigning the root of tree T E T<m to be of type 1 and the 
other vertices of the tree to be of type 2. Note that if M has maximum degree A, then every 
G E T<M also has maximum degree A. 

In this framework, Weitz's result establishing SSM for all graphs of maximum degree A 
when A < Ac(Ta) can be stated as establishing SSM with uniform rate for all G E T<m with 

M = ^ ^ ; and we are interested in establishing its analogy for general M. To this end, 

we will use the following notion of SSM for M. 

Remark 1. To establish SSM for M , it suffices to prove that SSM holds with uniform rate for 

all trees in T<m due to Corollary 2. In addition, note that SSM holds for M — if 

and only if it holds for (A) since the root of a tree T E J^<m is the only possible vertex of type 
1 in T . 

Finally, we define SSM for a branching matrix M . 

Definition 4. Given a branching matrix M , we say SSM holds for M if SSM holds with uniform 
rate for all G E J^<m- 
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Remark 2. To establish SSM for M , it suffices to prove that SSM holds with uniform rate for 
all trees in J^<m due to Corollary 2. 



3.2 Implications of SSM 

We present a new approach for proving SSM for a branching matrix M. There are multiple 
consequences of SSM for M as summarized in the following theorem. We first state some 
definitions needed for stating the theorem. 

Following Goldberg et al. [ ] we use the following variant of amenability for infinite graphs. 
Here we consider an infinite graph G — {V^E). For v and a non- negative integer d, let 
B^(i;) denote the set of vertices within distance < d from where distance is the length of the 
shortest path. For a set of vertices 5, the (outer) boundary and neighborhood amenability are 
defined, respectively, as: 

dS \— {w '. w ^ and w has a neighbor y ^ S} and r^ — sup ^!^^}} • 

The infinite graph is said to be neighborhood- amenable if inf^ r^ — 0. 

Now we can state the following theorem detailing the implications of SSM of interest to us. 

Theorem 3. For a t x t branching matrix M , if SSM holds for M then the following hold: 

1. For every G E J^<m, SSM holds on G. 

2. For every infinite graph G G T<m, there is a unique infinite- volume Gibbs measure on G. 

3. If M has maximum degree A, ift = 0(1) and A = 0(1), then for every (finite) G E J^<m, 
Weitz^s algorithm [ ] gives an FPAS for approximating the partition function Z(G). 

4. For every infinite H G J^<m which is neighborhood- amenable, for every finite subgraph 
G = (V^E) of H, the Glauber dynamics has 0{v?) mixing time. Moreover, if H — 
for constant d, then for every finite subgraph G — {V^ E) of H , the Glauber dynamics has 
O(nlogn) mixing time. 

Proof. Part 1 is by the definition of SSM for M. The uniqueness result follows from the fact 
that the infinite-volume extremal Gibbs measures on the infinite graph G can be obtained by 
taking limits of finite measures, see Georgii [14] for an introduction to infinite- volume Gibbs 
measures, and see Martinelli [ ] for Part 2. Part 3 immediately follows from the work of Weitz 
[3 ]. Finally, for Part 4, there is a long line of work showing that for the integer lattice in 
fixed dimensions, for the Ising model SSM on implies O(nlogn) mixing time of the Glauber 
dynamics on finite subregions of Z*^, e.ff., see Cesi [ ] and Martinelli [ ] (and the references 
therein) for recent results on this problem. These results for the Ising model are typically stated 
for a general class of models, but that class does not include models with hard constraints, such 
as the hard-core model studied here. Dyer et al. [lx] showed a simpler proof for the hard-core 
model that utilizes the monotonicity of the model. We use this result of [ ] in Theorem 8 
to get O(nlogn) mixing time for subregions of 1? . Goldberg et al. [ , Theorem 8] showed 
that for fc-colorings, if SSM holds for an infinite graph G that is neighborhood-amenable, the 
Glauber dynamics has O(n^) mixing time for all finite subgraphs of G. Their proof holds for 
the hard-core model which implies Part 4. □ 
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4 Establishing SSM for Branching Matrices 



In this section we present a sufficient condition implying SSM for the family of trees generated 
by a branching matrix. As a consequence of the approach presented in this section we get a 
simpler proof of Weitz's result [ ] implying SSM for all graphs with maximum degree A when 
A < Ac(Ta). We then apply the condition presented in this section to 1? in Section 5. 

To show the decay of influence of a boundary condition p, a common strategy is to prove 
some form of contraction for the 'one- step' iteration given in (1) below. More generally, we 
will prove such a contraction for an appropriate set of 'statistics' of the unoccupied marginal 
probability. 

A statistic of the univariate parameter x G [a, 6] is a monotone (i.e., strictly increasing or 
decreasing) function (p : [a, b] R. For di t x t branching matrix M we consider a set of t 
statistics (fi, . . . ,(ft, one for each type. For the simpler case when M = [A] and hence t = 1, 
we have a single statistic ip. Our aim is proving contraction for an appropriate set of statistics 
of the probability that the root of a tree is unoccupied. 

We first focus on the case of a single type. Consider a tree T = (1/, E) G T<m with root 
r. For V ^ V, let N{v) denote the children of v, and let d{v) := \N{v)\ the number of children. 
Let Ty denote the subtree rooted at v. We will analyze the unoccupied probability for a vertex 
V, but V will always be the root of its subtree. Hence, to simplify the notation, for a boundary 
condition p on. S CV, let — , ,. 

A straightforward recursive calculation with the partition function leads to the following 
relation: 

a^^{ ^ ''"^^'^^^ (1) 

^ I 1 I X -TT p otherwise. ^ ^ 

V l+^lUeiV(^)«u. 

Note, the unoccupied probability always lies in the interval / := 1 , i.e., for all all p. 

For V ^ let rriy := ^{a^) be the 'message' at vertex v. The messages satisfy the following 
recurrence: 

'^'^ ( 1 + A Y{^^N{v) «S ) ( 1 + A n^^^(^) ^-^{m^n 
Our aim is to prove uniform contraction of the messages on all trees T G J^<m- To this 
end, we will consider a more general set of messages. Namely, we consider messages mi, . . . , ttia 
where for every 1 < i < A, = and ai ^ I 1 . This set of tuples ai, . . . , a/\ G / 

contains all of the tuples obtainable on a tree. 

For ai, . . . , ^ let rui — ^{ai)^ 1 < ^ < A, and let 

/ 1 \ 

F(mi, . . . ,mA) := ^ 



.i + Anti^-n^.)y ' 

Ideally, we would like to establish the following contraction: there exists a < 7 < 1 such 
that for all , aA, Q^i, • • • , Q^a ^ ^' 

|F(mi, . . . , ttia) — F{mi^ . . . , rn/s^)\ < 7 max \mi — 

l<i<A 

where rui — ^{ai) and m[ — ^{oi[). We will instead show that the following weaker condition 
suffices. Namely, that the desired contraction holds for all \oii — Oi[\ < e for some e > 0. This is 
equivalent to the following condition. 
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Definition 5. Let I 



l+A' 



. For the branching matrix M = [A] , we say that Condition (^) 
is satisfied if for a// ai, . . . , a a ^ I, by setting mi = ^(o^i) for I < i < A, the following holds: 



\VF (mi, . . . , ttia) 



i=l 



dF (mi, . . . ,mA) 



drrij 



< 1. 



(^) 



Let us now consider a natural generalization of the above notion for a branching matrix 
with multiple types. Let M be a t x t branching matrix. For 1 < ^ < t, let = J2\=i 
denote the maximum number of children of a vertex of type £. Once again, consider a tree 
T = (y, E) G T<M with root r. For v ^ let t{v) denote its type. As before, N{v) are the 
children of d{v) is the number of children of v, and for a boundary condition p on. S CV^ aS 
is the unoccupied probability for v in the tree Ty under p. 

The recursive calculation in (1) for in terms of a^^w G N{v)^ still holds. For the case of 
multiple types, for v ^V^ let m^ := (ff(^y^{av) be the message at vertex v. The messages satisfy 
the following recurrence: 



mr 



For each type 1 < ^ < t, we consider contraction of messages derived from all o^i, . . . , a^^ ^ I- 
We need to identify the type of each these quantities ai in order to determine the appropriate 
statistic to apply. The assignment of types needs to be consistent with the branching matrix 
M. Hence, let 5^ : {1, . . . , A^} ^ {1, . . . , t} be the following assignment. Let <o = and for 
1 < i < let Mi^<i = J2l=i Mi,k- For 1 < i < t, for M^,<,_i <j< Mt^<i, let 5^(j) = i. 

For type 1 < ^ < t, for q;i, . . . , ^ set ruj — (fg^(^j^(aj)^ 1 < j < A^, and let 



F^(mi, . . . ,mAj := (fi 



i+Anfii^-;,)K) 



Note, 



, m; 



where N{v) = {^i, . . . , w^^v)}.^ 



(2) 



We generalize Condition (^) to branching matrices with multiple types by allowing a weight- 
ing of the types by parameters ci, . . . , q. 



Definition 6. Let I = 1 . For a t x t branching matrix M , we say that Condition (^^) 

is satisfied if there exist ci, . . . , such that for all 1 < i <t, for a// o^i, . . . , a^^ ^ I, by setting 
— ^si{i)(^i) for 1 < i < A^^ the following holds: 



A, 
1=1 



dFi (mi, . . . ,mA, 



The following lemma establishes a sufficient condition so that SSM holds for M. 



^Strictly speaking, requires arguments, so for (2) to hold in the case when d{v) < A^ we can simply 
add additional arguments corresponding to a = 1, which fixes these additional vertices to be unoccupied (and 
therefore absent). 
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Lemma 4. For a t x t branching matrix M , if for every 1 < £ < t, ipi is continuously dijfer- 



entiable on the interval I = 



l+A 



1 



and inf > 0^ and if Condition (i^) is satisfied for 



t — 1 or Condition fW^^ is satisfied for t > 2 then SSM holds for M , and hence the conclusions 
of Theorem 3 follow. 

Proof. For a tree T — iV^E) with root r, let and aj^^ denote the marginal probabilities 
that the root of T is unoccupied conditional on the vertices at level L (i.e., distance L from the 
root) being occupied and unoccupied, respectively. 

The main result for proving Lemma 4 is that there exist 7 < 1 and Lq < oc such that for 
every tree T G J^<m and every integer L > Lq, 



a 



a 



L,r 



< Y 



(3) 



We first explain why (3) implies Lemma 4 and then we prove (3). Consider a tree T = {V, E) 
with root r, and a boundary condition p on. S CV. Set L = dist(r, S) as the distance of S to the 
root of T. The hard-core model on bipartite graphs has a monotonicity of boundary conditions 
(c.f., [ ]) which implies that for odd ^ > ar > 

Hence, for any pair of boundary conditions p and r] on 5, 



1:^^, and for even L, a^^ > cxr > (^l^- 



< 



L,r L,r 



Therefore, by the definition of WSM in Definition 1, proving (3) implies WSM for T. Since 
this holds for all G J^<Mi by Observation 1, it imphes SSM for all V G J^<Mi which, by 
Remark 2, implies SSM for M. 

We now turn our attention to proving (3). Fix a, t x t branching matrix M and consider a 
tree T = (V^E) G T<m with root r. Given y G [0, 1], let ^L,v{y) denote the marginal proba- 
bility that the root of Ty is unoccupied given all of the vertices at level L (in Ty) are assigned 
marginal probability y of being unoccupied (conditional on its parent being unoccupied). Intu- 
itively, /3L,v(y) can be thought as the marginal probability conditioned on a 'fractional' boundary 
configuration at level L. As in (1), ^L,r{y) satisfies the following recurrence for y G [0, 1]: 



y 
1 

l+A 
1 



1+An 



Ay) 



if L = 0, 

if L > and N{r) 
otherwise. 



(4) 



/3L,r (1) and a 



From (4) and (1), it follows that aj^^ 
analyze the messages for a 
V let ruL^y (y) = y^t(v) {Pl,v (v))- Analogous to (2), we now have that: 



/3L,r(0). Hence, in order to 



and Q^^^, we will analyze the messages for ^L,r{y)- Therefore, for 



mL,r{y) = Ft^r) ( (y) , . . . , rriL- 



where N{r) = {^1, . . . , Wd(r)}' 



l+A 



1 



and hence we can 



Observe that for ah y G [0, 1], ah L > 0, ah G V, ^l,v (y) ^ I 
use Condition (^^) to analyze m^^r- 

Using the fact that PL,v{y) and niL^v (y) are continuously differentiable for y G [0,1], we 
have that for L > 0, 



a 



L,r 



a 



L,r 



= |/3L,r(l)-/3L,r(0)| < 



Jo 



dpL,r (y) 



dy 



dy < 



dmL,r{y) 



dy 



dy 



inf 
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By the hypothesis of Lemma 4, we know that ^[^^^ (x) > 0. Therefore, to prove the desired 
conclusion (3), it suffices to prove that there exist constants K < oo and 77 < 1 such that for 
every tree T G T<m with root r, ah L > 0, 



drriL^r {y) 



dy 



(5) 



Note that K and 77 should be independent of T and L, but may depend on A, (/^i, . . . , (/p^ and 
ci, . . . , Q. The constant K will be the following: 



K :-- 



max sup I (/^^ (^)l 
i<i<t xei 



mm cp 

i<i<t 



and the constant 77 will be the constant implicit in Condition (^). 

We will show (5) by induction on L. First we verify the base case L — 1. In this case, 



1 



Thus, 



druL^r (y) 



dy 



dy 



< sup 

xei 

< sup 



^t{r) ( 



X] 



sup 



< sup 



yem {1 + Xyd(r)f 
AA 



since L = 1 



by the chain rule 



< Ct(r)K 



by the definition of K. 



This completes the analysis of the base case. 

Now we proceed toward establishing the necessary induction step using the inductive hy- 
pothesis. We have that 



drriL^r (y) 



dy 



dFt^^) mL-i,wi (y) ^••-^ rriL-i^Wdfr) iv) 



dy 

^ dFt(r) {mi,.. . ,m^(^)) dmL-i,w^{y) 

i=l 



d(r) 
i=l 



drrii dy 
dFt{r) (^1, . . • , ^d(r)) 1 dmL-i,w,{y) 



where rrii := rriL-i^wiiy) 



drrij 



dy 
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drrij 



X max 



1 



l<i<d{r) 



dy 



by Holder's inequality. (6) 



From (^), there exists a universal constant 77 < 1 such that 



i=l 



drrii 



Therefore, it follows that 

drriL^r {y) 



dy 



< 77 q(^) • max 



1 



l<i<d{r) 
L-1 



dmL-i,w^{y) 



dy 



by (6) and the definition of 77 
by the inductive hypothesis. 



□ 



This completes the proof of (5), and hence that of Lemma 4. 
4.1 Reproving Weitz's Result of SSM for Trees 

In this section, we aim at finding a good choice of statistics. First we find such a statistic for 
the case iW = [A], i.e., the case of a single type, which enables us to reprove Weitz's result [ ] 
that when A < Ac(Ta) SSM holds for every tree of maximum degree A. 

Using Lemma 4 (and the simpler condition (^) for the case of a single type) we obtain a 
simpler proof of Weitz's result [ ] that for every tree T with maximum degree A + 1 (hence, 
for every graph G of maximum degree A + 1) and for all A < Ac(Ta+i) = A^/ (A — l)^^"*^, SSM 
holds on T (and on G). 

Theorem 5. Let (p{x) = ^ log (^^^^ where s = Then, Condition (i^) holds for M = [A] 

and A < Ac(Ta+i). Consequently, SSM and the conclusions of Theorem 3 hold for M — 

A + r 



A 



and A < Ac(Ta+i). 



Proof. First, a straightforward calculation implies that 



dF 



where ai — Lp ^(rrii) and a = (^1 + A fli^i 
Hence, we have 



dnii 

A 



1 — a 
s — a 

1 



(5 - ai), 



|vf||, = J: 

i=l 
A 



OF 



drrij 



1 — 

^-^ s — a 

1=1 
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< A I <5 ~ ( Y[ ) arithmetic-geometric mean inequality (7) 



We now use the following technical lemma. 
Lemma 6. 



(l-^)(l + s-(^i^y ^ 



max ^ — < 



where A is a positive integer and uo is the unique solution to uo{l + uo)^ — \. 
Using the above inequality (8) with Lemma 6, we have that: 

IIVFIL < 1 if • A < 1, 

where uj is the unique solution of uj{l + uj)^ = A. This leads to the desired condition A < 
Ac(Ta+i) = A^/(A - 1)^+1 so that SSM holds for M = [A]. As we noted in Remark 1, this is 

equivalent to SSM for iW = . This completes the proof of Theorem 5. □ 

Proof of Lemma 6. Let ^,,{x) = (^)i and f{x) = + Since = 

$A is a decreasing function in [0, 1] such that $a(0) = +oc and $a(1) = 0. Therefore it has a 
unique fixed point that can be shown to be x = Moreover, it is the case that $a(^) > a:^ if 

(l-\- — )(^(x)—x) 

and only if x < x. To prove Lemma 6, we notice that f{x) = 7^f/^T_^)2 ^ hence f{x) > for 
X < X and f^{x) < for x > x. This implies that / has a maximum at x, namely /(x) = □ 

4.2 DMS Condition: A Sufficient Criterion 

Theorem 5 suggests choosing (pj{x) = ^ log ^^-3^^ with appropriate parameters sj for a general 
branching matrix M. Under this choice, we obtain the following condition for SSM. 

Definition 7 (DMS Condition). Given atxt branching matrix M and A* > 0^ for 5i, . . . , 5^ > 1 

and c = (ci, . . . , > 0^ let D and S he the diagonal matrices defined as 

Djj — sup ^ — and Sjj = sj, 



where 



and A j = M^^ . 



We say the DMS Condition holds for M and A* if there exist si, . . . ,st > 1 and c > such that: 

(DMS) c<c. 
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Theorem 7. // the DMS Condition holds for M and A* > 0^ then Condition (i^) holds with 
the choice of ipj(x) = 7: log (^j^z^^ for all A < A*. Consequently, SSM and the conclusions of 
Theorem 3 hold for M and all A < A*. 



Proof. First, one can check that 



where ai — (f-^{mi) and a 



OF. 



druj 



Hence, it follows that 



1 — a 



Sn — a 



i=l 



dFj 



drun 



l — a 



Sj - a . 

1=1 



< 



< 



1 — a 



1 — a 



1 — a 



1 — a 



( A, 



1/AA 



i=\ 



\ 

( A, 

^^3^^3^- ^3 I n^^' 



1/A, 



\ — a 



1 — a\ 



1/A, 



1/A, 



\ — a 



by the arithmetic-geometric mean ineq. 

A 



l/A, 



by the definition of Qj 



^ji ^3i 



^ ^ A*a 

< Djj Mj£C£S£ by the definition of i^jj 

< Cj by the DMS condition. 

which satisfies the desired condition (^^) of Lemma 4. This completes the proof of Theorem 7. 

□ 



5 Application to in the hard-core model 

In this section, we show how to apply Theorem 7 and Theorem 3 to the two-dimensional integer 
lattice and improve the lower bound on Ac(Z^), resulting in the following theorem. 

Theorem 8. There exists atxt matrix M such that rsaw(^^) ^ ^<M cind the DMS Condition 
holds for A* = 2.3882. 

Therefore, the following hold for 1? for all A < A* ; 

1. SSM holds onZ^. 

2. There is a unique infinite- volume Gibbs measure on 1? . 
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3. If M has maximum degree A, if t — 0{1) and A = 0{1), then for every finite subgraph 
G of 1? , Weitz's algorithm [ ] gives an F PAS for approximating the partition function 
Z{G). 

4. For every finite subgraph G ofl?, the Glauber dynamics has O(nlogn) mixing time. 

We first illustrate our approach by showing that Theorem 8 holds with A* = 1.8801 for a 
simple choice of M. We then explain how to extend the approach to higher A. 

The graph 7? is translation-invariant, hence the tree Tsaw(^^5 ^) is identical for every vertex 
V ^1? . Fix a vertex, call it the origin o, and let us consider Tsaw(^^) = 7saw(^^5 0)- Each path 
from the root of TQQ^^^{7?) corresponds to a self-avoiding walk in 7? starting at the origin. Any 
walk on 1? starting at the origin o can be encoded as a string over the alphabet {N^ E, S, W} 
corresponding to North, East, South and West. The tree rsaw(^^) contains such strings, trun- 
cated the first time the corresponding walk completes a cycle. A relaxed notion of such a tree 
would be to truncate a walk only when a 4-cycle is completed. Denote such a tree by T4, and 
clearly we have that Tq^~s^{I?) is a subtree of T4. Our first idea is to define a branching matrix 
N so that T4 G J^<Ni and hence Tqq^^^{1?) G J^<n- 

To avoid cycles of length four, it is enough to track the last three steps of the walks. Labehng 
the paths using {N^ E, S, W} as mentioned above, their branching rule is easily determined. 
For example, a path labeled NWW is fohowed by paths labeled WWS, WWN and WWW 
which corresponds to adding the directions S, N and W respectively. As another example, a 
path labeled NWS is followed by paths labeled VKS'VK and WSS corresponding to adding the 
directions W and S to the path, while adding the direction E would have resulted in a cycle of 
length 4. The number of types in the corresponding branching matrix is < 4 + 4^ + 4^ < 5^. 
Indeed, we can reduce the representation of such paths by using isomorphisms between the 
generating rules among them. This results in 4 types in the following branching matrix N: 



N 



fo 


4 





o\ 





1 


2 








1 


1 


1 


{0 


1 


1 


0) 



(9) 



where the type z = 0, 3 of a vertex (walk) in the tree represents the fact that a continuation 
with a minimum of 4 — i additional edges are needed to complete a cycle of length 4. 

See Figure 2 for an illustration of this branching matrix N. One can verify that this branching 
matrix captures, inter alia, the self-avoiding walk trees from Z^: 

Observation 2. For any finite subgraph G = (V, E) of I? and v , Tqq^^{G^ v) G T<n- 

For this branching matrix, one can check that the (DMS) condition of Theorem 7 holds with 
A* = 1.8801, S = Diag(1.040, 1.388, 1.353, 1.255) and c = (0.266037, 0.100891, 0.100115, 0.0973861). 
Checking the DMS Condition for a given choice of parameters would have been a straightfor- 
ward task, were it not for the irrationality of the coefficients Djj. However, one can establish 
rigorous upper bounds for Djj, based on concavity of the function (of a) used in the definition 
of Djj^ in a suitable range of the parameters. These details will be discussed further below. As 
a consequence, we can conclude that Theorem 8 holds for N and A* = 1.8801. 

The primary reason why the branching matrix N improves beyond the tree-threshold of 
A < Ac(T4) = 27/16 = 1.6875 is that the average branching factor of any T G Tkn is significantly 
smaller than that of the regular tree of degree 4. 
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Figure 2: Assignment of the four types from matrix N defined in (9) to the self-avoiding walk 
tree Tgaw from Figure 1. In the circled area, we also draw redundant leaves at vertex j which 
may appear in the branching rule, but not in Tgaw 



To obtain a further reduction in the average branching, we observe that N did not consider 
the effect of occupying (or unoccupying) certain leaves as prescribed in Weitz's construction. 
Starting with T4, prune the leaves as is done in the construction of rsaw(^^) from Section 2.2. 
Denote the new tree as T!^. Clearly we still have that rsaw(^^) is a subtree of T4. 

Let us illustrate the difference between T4 and the pruned tree T4. We first fix an underlying 
order for the neighbors of each vertex. To this end, say N > E > S > W and this prescribes an 
ordering of the neighbors of each vertex. Consider a leaf vertex in the tree T4 corresponding 
to the vertex v in and to the path p in 7?. Since is a leaf vertex in T4, p must end with a 
cycle at say WNES. Since v was exited in the West direction at the beginning of the 4-cycle, 
and since W < the leaf vertex would be labeled occupied in Weitz's construction, thus 
resulting in the removal of and its parent in the construction of T4. Note, every vertex 
in T4 of type WNE has a child of type NES, and consequently (and its subtree) will be 
removed from the tree in the pruning process to construct T^. Thus, after removing vertices of 
type WNE (and similarly, WSE, SEN and ENW) from T4, it is stih the case that Tsaw(Z^) is 
a subtree of the resulting tree (T4). This highlights why T4 has a significantly smaller average 
branching factor than T4. 

We can define a branching matrix with 17 types (as illustrated in Figure 3), such that 
^4 ^ ^<M2i hence Tsaw(^^) ^ ^<M2- We can prove the DMS Condition is satisfied for M2 
at A* = 2.1625, as we will describe shortly, which significantly improves upon our initial bound 
resulting from considering T4. 

A natural direction for improved results is to consider branching matrices corresponding to 
avoidance of larger cycles, while also accounting for the removal of vertices prescribed by the 
construction of Weitz. We briefly outline such an approach for walks avoiding cycles of length 
at most 4, 6, and 8, respectively. Avoiding cycles of length 2i results in X]j<2i-i ^"^ — 
types, hence the computations become increasingly difficult for larger i. For 8-cycles the task of 
finding appropriate parameters to satisfy the conditions of Theorem 7 is still feasible. 

More precisely, we can define branching matrices Mi for i > 2, that (i) represent the structure 
of trees of walks avoiding cycles of length < 2z, as well as (ii) account for the removal of vertices 
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Typel Type 2 Type 3 Type 4 Type 5 Type 6 Type? Type 8 Type 9 



f I I ► ^ 1 I I ^ 1 I 1 

Type 1 Type 1 1 Type 1 2 Type 1 3 Type 1 4 Type 1 5 Type 1 6 Type 1 7 



Figure 3: Shapes that the seventeen types (or labels) represent for M2 where G J^<M2- 



based on children being labeled 'occupied.' One can extend the above construction of M2 for 
general i > 2 by using types encoded by longer paths with length at most 2i and ruling out the 
types that either contain a cycle of length at most 2i or whose children end up being labeled 
occupied. We can make the following general observation from our construction. 

Observation 3. For any finite subgraph G = iV^E) of I? and v , Tq^^{G^v) G T<Mi for 
any i > 2. 

As mentioned earlier, the matrix M2 constructed above consists of 17 types. An explicit 
description of it is shown in the Online Appendix [ )5], along with the associated parameters S 
and c for which one can check the DMS Condition for A* = 2.1625; this establishes Theorem 8 
for M2 and A* = 2.1625. 

The following table summarizes the threshold A* we obtain for each Mi'. 



Max length of Avoiding-cycles 


Effect of Occupations 


Number of Types 


A* 


4 


No 


4 


1.8801 


4 


Yes 


17 (< 


2.1625 


6 


Yes 


132 (< 5^) 


2.3335 


8 


Yes 


922 (< 5^) 


2.3882 



Note that, one can further improve the bound on A by using more types for higher i and 
hence Theorem 8 on will hold with the corresponding activity A*. For any such matrix, the 
verification of the DMS Condition relies on (i) 'guessing' appropriate values for the parameters S 
and c and (ii) formally verifying that DMS Condition holds for the chosen S and c. In choosing 
desirable S and c, we employed a heuristic random walk algorithm. 

To verify that the DMS Condition holds for a given rational matrix S and vector c is 
straightforward, provided we can obtain a rational upper bound for each type j for the function: 




Indeed, due to the concavity of this function for < 9j < 1, sj > 51/50 and A > 27/16, ^ it is 
always possible to find a provable upper bound for fj in such a regime. This can be done, for 

^This is a nontrivial algebraic fact. It can be proved by transforming the second derivatives condition to a set 
of integer polynomial constraints and using the "resolve" function in MATHEMATICA for the satisfiability of 
the constraints, which is rigorous by the Tarski-Seidenberg Theorem ] for the real polynomial systems [ ] and 
the so-called cylindrical algebraic decomposition [2]. 
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example, by describing a suitable 'envelope' for fj consisting of a piecewise linear function of 
the form: 

{Bi if a < a£ 

mm{bi {a - ai) + Bi, bu {a - au) + Bu} ii ai < a < au 
Bu ii a> au 

where a^^au are points such that b^ > fj{a£) > 0, bu < fj((^u) < 0, > fj(a£) and 
Bu > fj (o^u)- It is clear for any such function that gj(a) > fj{a), thus we obtain a provable 
upper bound for fj using gj. 

For every Mi in the above table, we provide S and c, along with appropriate envelopes that 
lead to upper bounds Djj for the corresponding Djj. Then we verify that the DMS Condition 
holds for the given values of A by replacing Djj with Djj. For z = 2, 3, 4 these values (M", 5', c, 
a£ and au) are given in the Online Appendix [35]. 



6 Ising Model 

The approach taken here for the hard-core model can also be employed to address corresponding 
questions in the well-studied Ising model. The Ising model, with inverse temperature parameter 
/3 > 0, on a finite graph G = (V, E) is the model associated with the Gibbs distribution fi on 
Vt = {-1, such that for a = [ai] E Vt, 

fi{a) = ^exp(^/3 aiaj^ , 

where the normalizing constant is the partition function: Z = Z{G, (3) Yla^VL ^^P (/^ j)eE 
The notions of SSM, the self-avoiding walk tree representation, and branching trees defined for 
the hard-core model extend identically to the Ising model (or, for that matter, any other 2-spin 
model). Moreover, an analog of Lemma 4 also follows easily for Ising. Then, by the use of an 
appropriate statistic (/p, the following simpler analog of the DMS Condition can be proved for 
the Ising model. 

Theorem 9. Given atxt branching matrix M and /3* > 0^ suppose there exists c— (ci, . . . , > 
such that 

tanh(/3*)Mc < c, (10) 
then SSM and the conclusions of Theorem 3 hold for M and all f3 G [0, ^5*]. 

Proof. First we note that Theorem 1 holds in general for all two spin models including the Ising 
model. Hence, Corollary 2 and Remark 1 are applicable to the Ising model as well. Further, 
observe that the proof of Theorem 3 (i.e., the implications of SSM) still hold for the Ising model. 
Consequently, we can prove Theorem 9 using similar notation and proof approach as was used 
for Theorem 7. 

Given a tree T G T<m and configuration p, let us define again a = ar{T^ j3) as the probability 
that the root r of T is minus-spinned. (Recall that in the hard-core model this was the probability 
that r was unoccupied.) li wi^ . . . are the children of r and Ti^ . . . ^T^ are the corresponding 
subtrees subtended at them, we let ai := aS^.(T^,/3) for i < k. For z > fc, we define ai := 
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1/2. Further let 



l-ai 



and I? 



l-a 
a 



Using these notations, a straightforward recursion 



calculation with the partition function leads to the following: 

A, 



L i}i + exp(2/3) 



(11) 



where j is the type of r, and Aj = Mj£. 

Motivated by (11), the function Fj (defined in Section 4 for the hard-core model) can be 
redefined for the Ising model as follows. 



Fj{mi, . . . ,mA. 



exp(2/3)(/^^./(mi) + l 



\ (/:^^./(m,) + exp(2/3) 



where ji is the type of child Wi and (fj is the statistic for a vertex of type j. Further, we define 
m := ^j('i9) and := It follows from (11) that m = Fj(mi, . . . ,mAj)- Then, one can 

prove the 'Ising version' of Lemma 4 with the interval / = [exp(— 2/3A), exp(2/3A)] using the 
same arguments as those in the proof of Lemma 4. Further, using the same arguments as in the 
proof of Theorem 7, with the choice of (pj{x) := log(x), we have that 



(e^^ - 1) 



< tanh {(3) , 



drrii {e'^Hi + 1) [e'^^ + i9i 
from which the desired condition (10) foUows easily. This completes the proof of Theorem 9. 

□ 

Using Theorem 9 with branching matrices M analogous to those we employed in Section 5 
for the hard-core model, we can prove that SSM holds for the Ising model on for all /3 < /3* 
as detailed in the following table: 



Dimension 


/?* 


2 


0.392190 


3 


0.214247 


4 


0.148045 


5 


0.113347 



In comparison, applying Weitz's general technique to I? implies SSM for (3 < .34657. 

We do not investigate the Ising model further because there are much stronger results known 
for this model. Onsager [ ] established that /3c(^^) = log(l + V^) ^ 0.440686. And for general 
trees, Lyons [ , Theorem 2.1] established the critical point for uniqueness. 
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